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Abstract 

A Bethe-Peierls treatment to dilution in frustrated magnets and spin liquids is given. 
A spin glass phase is present at low temperatures and close to the percolation point 
as soon as frustration takes a finite value in the dilute magnet model; the spin glass 
phase is reentrant inside the ferromagnetic phase. An extension of the model is given, 
in which the spin glass / ferromagnet phase boundary is shown not to reenter inside the 
ferromagnetic phase asymptotically close to the tricritical point whereas it has a turning 
point at lower temperatures. We conjecture similar phase diagrams to exist in finite 
dimensional models not constraint by a Nishimori's line. We increase frustration to study 
the effect of dilution in a spin liquid state. This provides a "minimal" ordering by disorder 
from an Ising paramagnet to an Ising spin glass. 
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1 Introduction 



When non magnetic sites are diluted in an unfrustrated ferromagnet with a probabiUty fi, the transi- 
tion temperature is reduced and vanishes at the percolation threshold 1 — /ip^o- In such non frustrated 
systems, two phases only exist: a low temperature ferromagnetic phase above the percolation thresh- 
old, and a paramagnetic phase. If the temperature is decreased at the percolation threshold, the dy- 
namics becomes slower because large-scale droplet-like objects of size form, with In^T ^ J/T 
These objects have energy barriers scaling like the logarithm of their volume ^. This results 
in a slow dynamics and interrupted aging Q (i.e. with a finite relaxation time |5|). This shows 
that despite the absence of frustration, the simplest models of dilute magnets already have a phe- 
nomenology close to the one of spin glasses, even though freezing in these systems is a cross-over due 
to an increasing correlation length becoming of order of the system size. This indicates that some 
perturbations of these unfrustrated systems may drive them to a true spin glass phase, which we 
show in the present article by studying the thermodynamics of a particular model. 

In dilute magnet compounds, such as Eua;Sri_a;S ||6|, |^, a low temperature spin glass phase 
appears close to the percolation threshold. The main features of the phase diagram are: 

(i) As the dilution ^ is increased from the pure system with = 0, the ordering temperature 
decreases. 

(ii) A tricritical point exists at a dilution and temperature T^, with \ — of order of the 
percolation threshold 1 — /ip^o in the absence of frustration. At this tricritical point, the 
ferromagnetic, paramagnetic and spin glass phases meet. 

(iii) As dilution is increased from /i^, the spin glass transition temperature decreases from Tt at 
the tricritical point to zero at the percolation threshold 1 — fip, with 1 — //p the percolation 
threshold of the system with frustration, smaller than the percolation threshold 1 — fip^o in the 
absence of frustration. 

(iv) The spin glass phase is reentrant inside the ferromagnetic phase. 

One purpose of the present article is to show that these qualitative features of the phase diagram can 
be reproduced in a model that combines dilution and short range frustration. This model does not 
consist in a detailed microscopic modeling of Eua;Sri_a;S, but rather contains the generic ingredients 
entering the physics of these systems (dilution and short range frustration) Q. Even though this 
treatment relies on a specific lattice topology (a tree structure), Bethe-Peierls phase diagrams are 
equivalent to mean field phase diagrams while the Bethe-Peierls method is powerful enough to give 
an exact answer to the issue of reentrance. The resulting phase diagrams are therefore not expected 
to be specific to our treatment but are generic features of the coupling Hamiltonian. 
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The article is organized as follows. The model is given in section g The paramagnetic phase 
boundary is solved in section ^. We study in section ^ the spin glass / ferromagnet phase boundary 



and show that it is reentrant. The issue of reentrance is non trivial because Nishimori's argument |10 



11, 12, 13, 14 1 does not hold in our model, part of the exchanges being frozen. An extension of the 
model is given in which the spin glass / ferromagnet phase boundary is shown not to reenter in the 
ferromagnetic phase asymptotically close to the tricritical point, whereas it has a turning point at a 
lower temperature. We expect this unusual type of phase diagram to be a generic feature of models 
combining disordered and frozen exchanges, and may be obtained in finite dimensional models also. 
Finally, we study in section |5| the effect of diluting a spin liquid state, in which case ordering by 
disorder generates a transition from an Ising paramagnet to an Ising spin glass. This ordering by 
disorder mechanism is "minimal" because spins have an Ising symmetry only. 

2 The model 

2.1 The dilute magnet model and its generalization 

In a Bethe-Peierls calculation, only the properties of the "top" spin (the highest one in the hierarchy; 
see Fig. |l|) are considered, and the thermodynamic limit is obtained by growing the number of 
generations to infinity []T5| , 16, The top spin fixed point magnetization distributions P*{m) are 



of three types: (i) paramagnetic phase: P*{m) = 5{m); (ii) spin glass phase: P*[m) is even; (iii) 
ferromagnetic phase: P*[m) has a finite first moment. The transitions between the phases are of a 
mean field type [15, |16|, |l^. The phase diagram of the ±J model on the Bethe lattice |l^, 20, 21 



22, 23, p3] is very similar to the one of the Sherrington-Kirkpatrick model [?], and the Bethe-Peierls 



treatment allows a correct description [gj] of the Almeida-Thouless line |25]. This shows that a 
Bethe-Peierls treatment succeeds in reproducing the mean field phase diagram of spin glass models. 

We consider a model in which the ferromagnetic bonds J of the Cayley tree are canceled with a 
probability ^. Frustration is added by completing the triangles, with frozen antiferromagnetic bonds 
r forming a Husimi cactus-like structure |26| (see Fig. ||). We use the binary variables = 0, 1, with 



Ji = OiJ. The temperature is expressed in units of J. The Hamiltonian is 

where {i,j) denotes the bonds of the tree structure and {i,j)' the next nearest neighbor pairs of sites 
in the same generation. The distribution of the ^-variable is 

p[e) = {l-^,)5{e-l) + ^i6{e). (2) 

The specificity that some bonds are frozen in this model has drastic consequences on the shape of 
the phase diagram, as we will show. 



2 



ferromagnetic coupling (9=1) 

— — — cancelled ferromagnetic coupling (9 = 0) 
antiferromagnetic coupling ( x ) 

Figure 1: The Husimi cactus- like structure model of dilute magnet with frustration. The structure 
with a top spin x is obtained by gluing the two structures with top spins y and z. The top spin is 
the highest one in the hierarchy (at site x on this figure). The ferromagnetic bonds of the tree are 
canceled with a probability fi and a fixed antiferromagnetic coupling r is added. 



The Husimi cactus structure allows the introduction of a local frustration resulting from next- 
nearest-neighbor interactions, and can be used to mimic the effects of local antiferromagnetic in- 
teractions in dilute compounds. Chandra and Dougot considered a frustrated spin model on 
a regular Husimi cactus structure, and studied ordering by disorder in the spin liquid state in the 
Bethe-Peierls limit (see also ||2^ for a study of the effect of quantum fluctuations). These authors 
considered a non disordered model in which a spin glass phase cannot exist |27|. In our model with 



randomness, we show the stability of a spin glass solution in some regions of the phase diagram. 

We consider bond instead of site percolation because the site percolation threshold of the Husimi 
cactus structure would be equal to the one of the tree structure, independent of the additional bonds 
r. The bond percolation model is therefore better suited for modeling dilute compounds |^], since the 
bond percolation threshold of the structure without frustration (a tree structure) is 1 — fj-p^ = 1/2, 
larger than the bond percolation threshold 1 — fip = 1 — l/\/2 of the structure with frustration (the 
Husimi cactus structure shown on Fig. ||). 

For the sake of generality, we not only consider a dilute magnet model with 6 = 0,1, but extend 
the bond distribution (^) to incorporate possible antiferromagnetic bonds on the tree structure. The 
distribution of the bond variables 6 is 

p{e) = (1 - A)(i - n)5{e - 1) + ii6{e) + a(i - /i)<5(0 + 1), (3) 

while the additional antiferromagnetic bonds r are frozen. This model interpolates between the dilute 
model with a short-range frustration r (A = 0), and the it J model (// = and r = 0). 
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2.2 Absence of a Nishimori line argument 

Our Hamiltonian is formally invariant under local gauge transformations ai eiai, Ji 



with ej = ±1 |29]. In some spin glass models (such as the it J model), gauge invariance provides strong 
constraints on the phase diagram. The internal energy can be calculated exactly on Nishimori's line 



by expanding the average energy over the gauge group [|10|. This line crosses the phase boundary 



at the tricritical point |10, 11, 12, Moreover, spin correlations can be related to gauge variable 



correlations, with the consequence that the frontier between the ferromagnetic and spin glass phases 



is either vertical or reentrant |1C]. In our model, gauge invariance is useless for the following reason. 
One can define a local distribution of bond variables Pij{J), being (|3|) on the tree bonds, and 
6{J + t) on the antiferromagnetic bonds r. Nishimori's line is defined by Pn = ^In^^^^^ [14|, 
and /3iv = +oo. The first equality originates from the tree bond variables and the second from the 
frozen antiferromagnetic bonds r. The two equalities can be formally met if A = 0, in which case 
Nishimori line is = +oo. However, this does not make predictions on the phase diagram possible 
even in the case A = [^] because one does not expect to be able to describe finite temperature spin 
glass properties in terms of the ground state only (the only state selected if = +oo). Nishimori's 
argument can therefore not be made in this model, and the question of reentrance (item (iv) in 
the introductory section) cannot be answered on the basis of Nishimori's line while Bethe-Peierls 
calculations are powerful enough to allow the derivation of exact results. We show that the spin 
glass / paramagnet boundary is reentrant in the dilution model (A = 0) in the {fj,, T) plane. In the 
=bJ model with the additional coupling r (^u = 0), and in the (A,T) plane, we show that the spin 
glass / ferromagnet phase boundary is not reentrant asymptotically close to the tricritical point, 
whereas it has a turning point at lower temperatures. This behavior, richer than in usual spin glass 
models, is to our opinion a generic feature of Hamiltonians combining disorder and frozen bonds. We 
conjecture the existence of finite dimensional models with a similar phase diagram. 



3 Recursion relations 

We now derive the recursion of the top-site magnetization when cacti are glued as shown on Fig. |l|. 
We denote by the magnetization at site x, and ruy and ruz the magnetizations at the descendant 
sites y and z. The derivation of the recursion relations in our dilute magnet model is similar to the 
case of the ±J model (see Ref. p^). 

3.1 Recursions and the paramagnet phase boundary 

Following Ref. |^3|, we denote by z'^^ the conditional partition function with the spin at site x 
frozen in the direction it. The magnetization at site x is nix = 

partition functions Z^'' are related to the partition functions z'^) of the descendant sites according 
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to 

Z^) = J2 Wls^{a,\ay,az)Z^'^y^Zt\ (4) 



with the Boltzmann weight factor 

W^y,eS^x\^y^ = ^^P iPi^yf^xf^y + dzCTxC^z)) exp {-Praya^), (5) 

and 6 = 0, ±1. We next trace over the spins at sites y and z in Eq. ^ to obtain 

m /) X my{9y -uOz) +mz{Oz -uOy) 
mx = f{my,mz\ey,ez) =p- L n ^ 1 2n n ^ ' (6) 

with p = tanh (/3J) and u = tanh (/?t). The recursion of the magnetization distribution is 

Pn+iimx) = J drnxdniy ^ p{9y)p{dz)Pn{my)Pn{mz)5 {nix - f{my,mz\6y,9z)) , (7) 

9y,9z 

with p(6) the distribution of bond variables and P„ the magnetization distribution of the top 
spin with n levels of hierarchy. We denote by {{m^))n the moment of order k of Pn{m). 

We now parametrize the tricritical line, where the three phases (paramagnetic, ferromagnetic 
and spin glass) meet. The meeting point of these phases is a line in the parameter space (A,^, T). 
If A [fi] is fixed and the phase diagram is considered in the (/i, T) [(A, T)] plane, the three phases 
meet in a tricritical point. Let us first consider the stability of the paramagnetic solution with 
respect to perturbations in the first moment. To lowest order the recursion of the first moment 
is {{m))n+i = 2p{{Gy^z)){{'m))n, with Gy^z = {Oy — 'u6z)/{l — up^OyOz)- The disorder average of 
G is understood as {{Gy^z)) = J2ey,ezPi^y)Pi^z)Gy^z- The paramagnetic solution is stable with 
respect to perturbations in the first moment if 2p{{Gy^z)) < 1- A similar reasoning shows that the 
paramagnetic solution is stable with respect to perturbations in the second moment if 2p'^{{Gy^z)) < 1- 
To summarize, the tricritical line is defined by 

2p{{Gy,z)) = 1 , and 2p^{{GlJ) = 1. (8) 

3.2 Limiting cases 

The pure system: Let us consider the recursion (I) in the pure system limit in which the variables 9 
are all equal to unity. This amounts to specializing the distribution (^) to the case A = ^ = 1 while 
keeping finite the local frustration r. The only possible phases are ferromagnetic and paramagnetic. 
The recursion of the magnetization is m„+i = 2p{l — u)mn / {1 — pu"^ + {p"^ — u)m'^) . The paramagnetic 
phase is stable against ferromagnetic fluctuations if 2p{l — u)/{l — pv?) < 1. The phase diagram is 
shown on Fig. |2| as a function of the frustration r with a spin liquid phase if r > 1. When considering 
in the following the frustrated magnet model, we assume r < 1, in which case the pure system has an 
ordered phase at low temperature. Diluting the spin liquid state with r > 1 is examined in section |5|. 
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Figure 2: Phase diagram of the pure Husimi cactus system (all the coupling 6 being unity) as a 
function of the frustration r. A low-temperature ferromagnetic phase is present if r < 1. If r > 1, 
the system is a spin liquid (it does not order even at zero temperature). 
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Figure 3: Phase diagram of the dilute magnet model with frustration (A = 0, r = 0.1). The spin 
glass phase exists below the percolating dilution fip = 2~^/^ ~ 0.707. The paramagnetic / spin glass 
phase boundary inside the ferromagnetic phase is unphysical, as well as the paramagnetic / spin glass 



boundary inside the spin glass phase. The solid line is obtained from the calculation in section 
of the frontier between the spin glass and ferromagnetic phases. This solution is exact close to the 
critical point and we have continued it to lower temperatures by an arbitrary linear behavior. The 
exact zero temperature spin glass / ferromagnet phase boundary is no ~ 0.24191. 
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lb J model: The it J model is recovered if = r = 0, with a tricritical point at coordinates pt = 1/2, 

Dilute magnet with frustration: If A = and r 7^ the relations become 

2p{l - - u){l - fiup^) ^ ^ 
1 — up"^ 

2(1 . 



^1 — up' 



2\2 



(1 - ^)(1 - ^/)^ + ^(1 + n^)(l - up'Y = 1, (10) 



determining the paramagnetic phase boundary shown on Fig. ^. The frontier between the spin glass 
and ferromagnetic phases will be examined in Section ^ by looking for an instability in the first 
moment of the spin glass solution. 

The zero temperature limit of the paramagnetic/spin glass phase boundary can be obtained by 
considering first the limit p = 1 in Eq. (|10|), and second the limit u = 1. The order in which the two 
limits are taken is imposed by the fact that 1— p<Cl — M<Clat low temperatures because r < 1. 
One finds the limit tobel — /i = l — 1/ \/2. As it is expected, this dilution is equal to the percolation 
threshold of the Husimi cactus structure. 

4 Frontier between the spin glass and ferromagnetic phases 
4.1 Method 

In order to determine the frontier between the spin glass and ferromagnetic phases, we study the 
instability of the spin glass solution with respect to perturbations in the first moment. This involves 
first calculating the spin glass solution close to the tricritical line and next determining whether 
this solution is stable with respect to ferromagnetic fluctuations. Carlson et al. |^3[ performed 
this calculation for the itJ model close to the tricritical point, and shown the spin glass phase 
to be marginally reentrant inside the ferromagnetic phase. By marginal, we mean that the spin 
glass / ferromagnet phase boundary has a quadratic behavior Af — A ~ (T^ — T)^, which is specific to 
this model. Other models (as the one we presently analyze) have a linear behavior Aj — A ~ — T, 
with a positive (reentrant behavior) or negative prefactor (non reentrant behavior). The calculation 
follows Ref. [^] where the it J model was solved, and is asymptotically exact close to the tricritical 
point. This is complemented by an exact determination of the zero temperature phase boundaries 
which, to our knowledge, has not appeared previously in the literature even for the it J model. 

We first determine the asymptotic spin glass solution close to the tricritical line. The second 
moment to lowest order is 

«" » = 4p?((0.,.0„„g». - 
with H = (p^Oye, - u)/(l — up OyOz), and the subscript "t" denoting a quantity evaluated on the 
tricritical line. We next consider a perturbation in the first moment of the spin glass solution. The 
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Figure 4: Possible shapes of the spin glass / ferromagnet phase boundary in the ±J model with 
a small antiferromagnetic coupling r corresponding to a spin glass phase not reentrant inside the 
ferromagnetic phase asymptotically close to the tricritical point. This implies two possible behaviors: 
(a): no reentrance at any temperature (b): no reentrance close to the tricritical point but reentrance 
at lower temperatures. We prove that (b) is the correct behavior in a zero temperature exact solution. 

recursion of the first moment is {{m))n+i = K{{m))n , with k = 2p{{Gy^z)) — '^Pt{{'ni'^)){{Gy^zH))t 
to order (Tf — T). If k < 1 the spin glass phase is stable and otherwise the ferromagnetic phase is 
stable. 



4.2 ±J model with an additional short-range coupling r — Fixed (A, T) phase 
diagram 

In the small-r limit, it is straightforward to show that 

1. The tricritical point can be determined in an expansion in ut- pt = 1 — 2Xt = "^(1 + 

2. The slope of the spin glass / ferromagnet phase boundary at the tricritical point is 

Ut ( , 5\/2^ 



dT 



1 



-0.467 Ut <0, 



(12) 



with = 1/ tanh ^ (l/-v/2) — 1.135 the tricritical point temperature with r = 0. 



/.From what we deduce that the spin glass phase does not reenter inside the ferromagnetic phase 
close to the tricritical point. 

Notice that dK,/dT\^ in Eq. (12) vanishes if r = 0. This is because reentrance is marginal in the 



±J Bethe lattice spin glass [23| and can therefore not be obtained from an expansion to first order 
in Tt — T. We have evaluated numerically the coefficient n with a finite r and a finite /i. As ^ is 
increased above a critical value, the transition changes from reentrant to non reentrant. 

We now derive the exact spin glass solution in the zero temperature limit, which allows to discrim- 
inate rigorously between the two behaviors on Fig. ^ (a) and ^ (b). We look for the zero temperature 
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fixed point spin glass and ferromagnetic solutions P*{m) under the form 

= ^^5(m - 1) + (1 - x)6{m) + ^^<5(m + 1), (13) 
with X and y the spin glass and ferromagnet order parameters. It turns out that the functional form 



of the magnetization distribution (13) is stable when it is iterated in the zero temperature limit of 
(^) and To determine x and y, we impose ( |l3| ) to be the fixed point magnetization distribution. 
The solution with a finite magnetization is x = (1 - 4A)(1 - 2A), and = (1 - 4A)(1 - 8A)/(1 - 2A)^ . 
Imposing > leads to the intersection Aq = 1/8 = 0.125 of the spin glass / ferromagnet phase 
boundary and the zero temperature axis in the (A, T) plane. Aq is independent of the strength of the 
additional coupling r. If r = 0, the value A^ of A at the tricritical point is A^ = \{\ — ■^) ~ 0.146 [23|, 
larger than Aq- The zero temperature solution is therefore consistent with the reentrant behavior 
of the spin glass / ferromagnet phase boundary of the itJ Bethe lattice spin glass pf. As r IS 



increased, the tricritical point (At(r), (r)) evolves with r whereas the intersection of the spin 
glass / ferromagnet phase boundary and the zero temperature axis remains equal to Aq. Therefore, 
if T is small, \t{T) remains larger than Aq. From what we deduce the existence of a turning point 
in the spin glass / ferromagnet phase boundary (Fig. ^ (b)): the spin glass phase does not reenter 
close to the tricritical point whereas it reenters at lower temperatures. We believe this behavior to 
be generic of spin glass models with frozen exchanges and we conjecture that a similar behavior may 
be obtained in finite dimensional models. 

4.3 Dilution with a short-range frustration r — A = 0; (//, T) phase diagram 

A small-r perturbation calculation leads to pt = 1 — ^ut, = \ — ^ut, and dK/dT\^ = t /2Tf > 0, 
which proves that the spin glass phase reenters in the ferromagnetic phase close to the tricritical point 
in the limit of small r. We have shown on Fig. |3| the behavior of the spin glass / ferromagnet phase 
boundary. This phase boundary is exact only close to the tricritical point, and we have continued 
it by an arbitrary straight line at lower temperatures. The reentrant behavior is confirmed by zero 
temperature exact results. The paramagnetic / spin glass frontier intersects the zero temperature axis 
at the percolation threshold, and the spin glass / ferromagnet frontier intersects the zero temperature 
axis at f^o, the real root of — lO/Ug + 6^o ~ ^^o + 1 = 0, approximately /Uq — 0.24191. This confirms 
the reentrant behavior of the spin glass transition in the dilute magnet model with frustration. 



5 Diluting the spin liquid 

We now consider dilution in the regime r > 1, i.e. when the pure system is a spin liquid (see Fig. |2|). 
As it could be expected, a ferromagnetic instability of the spin liquid solution (Eq. |9|) does not exist. 
However, a spin glass instability of the paramagnetic solution does exist upon diluting the system. 
Let us first consider the zero temperature phases, in the limit 1 — u<^l— p<^l (since r > 1). The 
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Figure 5: Boundary between the spin glass and paramagnetic phases upon diluting the spin liquid 
state of our model (r > 1). The spin glass phase is confined inside the boundary shown for r = 1.2 
(O), r = 1.25 (+), and r = 1.3 (□). The spin glass boundary collapses onto the point // = 1/2 
in the zero temperature limit. This boundary behaves like T ~ l/ln|/u — 1/2| around this singular 
point The spin glass phase is favored upon increasing the temperature (ordering by disorder - see 



section 6.2) 



phase diagram at finite temperatures is shown on Fig. |5[ A finite temperature spin glass phase opens 
from the point = 1/2, T = 0) as temperature is increased from zero. The low temperature phase 
boundary is T = 4(1 — T)/[ln (^ — 1/2)^]. This provides a simple situation in which diluting an Ising 
spin liquid results in an Ising spin glass phase. The underlying ordering by disorder mechanism |36| 



is analyzed in section 6.2 



6 Conclusion: diluting a frustrated magnet versus diluting a spin 
liquid 

We have shown the existence of a spin glass solution upon diluting both the weakly frustrated magnet 
(r < 1), and the spin liquid (r > 1). We underline the differences in the physics in these two regimes. 

6.1 Diluting the frustrated ferromagnet 

We believe the generation of a spin glass phase upon weakly frustrating a dilute magnet close to the 
percolation threshold to be due to the following: the strong diluted unfrustrated magnet is already 
close to a spin glass. This can be seen on the example of square lattice dilute magnets [^], where 
dilution removes sites in the ferromagnet up to the point where the percolating cluster is a fractal 
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object at the percolation point. Since the order of ramification of percolating clusters is finite [32|, 
one can isolate large droplet-like objects from the remaining of the structure by cutting a finite 
number of bonds. This results in large sets of spins that can be reversed at a finite energy cost, thus 
being responsible for the existence of quasi-degenerate ground states separated by a large distance in 
phase space (with different magnetizations |^^), and with barriers scaling like the logarithm of their 
volume Q. 

The addition of frustration in dilute magnets close to the percolation threshold turns the quasi 
spin glass order into a true one. We have shown this explicitly in our model, and a similar behavior 



was obtained in another model [12|. We do not expect the low energy states of the spin glass phase 
with a small frustration to be very different from the droplet-like states of the unfrustrated magnet. 

The "chaos and memory" behavior of metallic spin glasses was put forward in Ref. [^], associated 
to the growth of fractal droplets with a chaotic behavior in the sense that droplets at a given 
temperature overlap weakly with droplets at a different temperature. We do not expect a chaotic 
behavior in our model because the finite temperature droplet excitations of the frustrated dilute 
magnet should be obtained from reversing clusters of spins in the unfrustrated magnet dilute lattice. 

6.2 Diluting the spin liquid: ordering by disorder 

The mechanism for generating a spin glass phase from the spin liquid is different. The spin glass 
phase originates from a balance between the small-dilution regime in which dilution suppresses the 
liquid behavior in favor of spin-glass correlations, and a large-dilution regime in which dilution 
suppresses spin glass correlations by cutting the system into finite pieces. This is an order by 



disorder mechanism [36|: thermal fluctuations favor a spin glass arrangement and therefore reduce 
the phase space dimensionality compared to the one of the spin liquid state. Let us think in terms of 
low temperature properties in the large-r limit. In this limit, the neighboring spins coupled by the 
strong antiferromagnetic exchange r correlate antiferromagnetic, thus leaving mainly two residual 
degrees of freedom per bond r. We note my and the magnetization of these two spins, and, for 
the sake of a qualitative argument, assume my = —mz as a result of the strong bond r. Let us 
assume the spins at sites y and z to be frozen and look whether freezing is relevant in the Bethe- 
Peierls limit. We see from Eq. (P) that = if the two ferromagnetic bonds are present (i.e. 
if the triangular plaquette is frustrated. By = 6^ = 1). In the unfrustrated plaquettes 9y = 0, 
0^ = 1, or = 1, ^2 = 0, correlations in the magnetization can propagate from one generation to 
the other. The system is cut into two pieces if ^j, = ^2 = 0, preventing correlations to propagate 
from one generation to the other. When the ferromagnetic bonds are diluted, frustration is reduced 
since the fraction of frustrated triangular plaquettes (1 — /i)^ decreases upon increasing the dilution 
ji. Decreasing frustration therefore decreases the short range liquid-like correlations and favors a 
cooperative spin glass arrangement. In the large dilution limit, the exchanges are severely depleted 
and a paramagnetic behavior is restored since the system is cut into finite pieces. In between these 
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two regimes, the unfrustred bond configurations By = l,9z = and 9y = 0, 9z = 1 with a weight 
2//(l — n) dominate the physics and make a spin glass order possible. 

6.3 Concluding remarks 

Finally, we would like to compare the present work to other approaches developed previously in the 
literature, and mention some open questions. The phase diagram with all the bonds drawn from the 



distribution (y) was studied by Aharony jS^] , Giri and Stephen and Viana and Bray ||3^ . These 



models share similarities with the dilute fee antiferromagnets studied by de Seze ||40[] and Wengel, 



Henley and Zippelius ||4l[| . 



Nieuwenhuizen and Nieuwenhuizen and van Duin studied the field theory of a model of site- 



disordered magnet 43 1. One may also define a model similar to ours in a finite dimension. As we 
conjectured, a phase diagram similar to the one on Fig. ^ (b) may be obtained. On the other hand, 
it may be useful to investigate replica symmetry breaking in Bethe-Peierls calculations. 

Hierarchical lattices have been used previously by Georges and Le Doussal ||l^ in relation with the 
renormalization group flow along Nishimori's line, and by Gingras and S0rensen to study reentrance 
from a paramagnetic to a ferromagnetic phase [^4|. A model with frozen exchanges may be studied on 
a finite dimensional hierarchical lattice, which could be a first step in addressing the phase diagram 
on Fig. 0(b) in a finite dimension. This approach should probably rely on a numerical iteration of 
the renormalization equations similar to Ref. |4^] while an analytic study was possible in the present 
work. 



Finally, ordering by disorder seems to be a generic behavior of spin liquids |2^, 28, 45, 46 1. 
We found in the present work an ordering by disorder resulting in a transition from a paramagnetic 
to a spin glass ordering in an Ising model. This may be viewed as a "minimal" ordering by disorder 
from a Z2-symmetric paramagnet to a spin glass because the Ising order parameter has the lowest 
possible spin symmetry. 
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